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We present a unified transport theory of hybrid structures, in which a confined normal state 
(N) sample is sandwiched between two leads each of which can be either a ferromagnet (F) or a 
superconductor (S) via tunnel barriers. By introducing a four-dimensional Nambu-spinor space, 
a general current formula is derived within the Keldysh nonequilibrium Green function formalism, 
which can be applied to various kinds of hybrid mesoscopic systems with strong correlations even in 
the nonequilibrium situation. Such a formula is gauge invariant. We also demonstrate analytically 
for some quantities, such as the difference between chemical potentials, superconductor order pa- 
rameter phases and ferromagnetic magnetization orientations, that only their relative value appears 
explicitly in the current expression. When applied to specific structures, the formula becomes of the 
Meir-Wingreen-type favoring strong correlation effects, and reduces to the Landauer-Biittiker-type 
in noninteracting systems such as the double-barrier resonant structures, which we study in detail 
beyond the wide-band approximation. We find that the spin-dependent density of states of the 
ferromagnetic lead(s) is reflected in the resonant peak and resonant shoulder structure of the FV 
characteristics of F/I/N/I/F structures with large level spacing. Fhe tunnel magnetoresistance 
that exhibits complex behaviors as a function of the bias voltage, can be either positive or negative, 
suppressed or enhanced within the resonant peak region(s), depending on the couplings to the leads. 
The Andreev current spectrum of F/I/N/I /S structures consists of a series of resonant peaks as 
a function of the gate voltage, of which the number and amplitude are strongly dependent on the 
bias voltage, degree of spin polarization of the ferromagnetic lead, energy gap of the superconduct- 
ing lead, and the level configuration of the central region. In S/I /N/I/S resonant structures with 
asymmetric superconducting energy gaps, the Josephson current through a single resonant level is 
slightly enhanced in contrast to the significant enhancement of the Josephson current in S/N/S 
Junctions. The current-phase relation is relevant to the level position and the couplings to the 
superconducting leads. 

PACS numbers: 72.10.Bg,72.25.-b,74.50+r,73.21.-b 
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L INTRODUCTION 

Electronic transport in mesoscopic systems or 
nanoscale structures has received extensive theoretical 
and experimental attention^ In mesoscopic systems the 
sample size is smaller than the phase coherent length, 
and electrons retain their phase when travelling through 
the sample. In the ballistic limit, i.e., when the dimen- 
sions of the sample are smaller than the mean free path, 
electrons can traverse the system without any scatter- 
ing. In contrast to macroscopic systems, the conductance 
of mesoscopic systems is sample-specific, since electron 
wavefunctions are strongly dependent on the form of the 
boundary of the sample and the configuration of scatter- 
ers located within the sample. 

To calculate the conductance of mesoscopic systems, 
one should first consider the wave nature of electrons. 
The classical Boltzmann transport equation^ is obvi- 
ously inappropriate, since the assumption that electrons 
can be viewed as classical particles does not hold at 
a mesoscopic scale due to the Heinsenberg uncertainty 
limitation. Linear-response theory^ is restricted to the 
weak perturbation regime and apparently can not be 



applied to the nonlinear or nonequilibrium situation. 
Electronic transport through a mesoscopic medium is in 
effect a wave transmitting process of electrons, which 
can be associated with a scattering matrix. In mea- 
suring the conductance one always connects the sam- 
ple to electron reservoirs through perfect leads^ In a 
two-terminal setup (£=lcft, 7?.=right), the Landauer- 
Biittiker formula^ states that the current X can be ex- 
pressed as a convolution of the transmission probability 
T and the Fermi distribution function f a (a = £, 1Z), i.e., 
1= 22/ T(e)[fc(e)-f n (e)]de. The conductance Q in the 

linear-response regime is Q = 2f_ J T(e)(—^)de. Such 
a formulation seems more appealing since the transport 
properties are encoded in the corresponding transmission 
probability, which can be calculated by various methods. 

The nonequilibrium Green function (NEGF) 
approachSiLSiS has proven to be a powerful technique 
to investigate transport problems in the many-body 
systems and mesoscopic systems. The equation of 
motion for the NEGF G < , the quantum Boltzmann 
equation (QBE)^ serves a starting point for many 
transport calculations in the many-body problems^^ 
where a four-variable distribution function is required to 
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incorporate the quantum effect due to the uncertainty 
principle. The Keldysh formalism of the NEGF)2ii£ 
due to its integral form, becomes a popular method in 
the formulation, calculation and simulation of recent 
mesoscopic transport problems. Caroli at al. were the 
first to employ the Keldysh NEGF technique to study 
the tunneling problems of a biased (nonequilibrium) 
metal-insulator-metal junctional Meir and Wingreeni^ 
in 1992 derived a useful formula for the current through 
an interacting region with normal leads and applied it 
to investigate the transport properties of a quantum dot 
in Kondo and fractional quantum Hall regimes. Later 
the Keldysh NEGF formalism was used to analyze the 
I-V characteristics of superconductor-superconductor 
point contacts and the transport problem in a quantum 
dot with superconductor leads in Nambu space— ii^i 5 - 
By introducing Green functions in the spinor space, the 
Keldysh NEGF approach has been also employed to 
study a quantum dot connected to two ferromagnetic 
electrodes— Therefore, incorporating both Nambu 
and spinor spaces is a convenient device in order 
to investigate transport problems in the presence of 
both superconductors and ferromagnets within the 
Keldysh NEGF formalism. It is the purpose of this 
work to present a unified theory of electronic transport 
through an interacting region connected to either bulk 
ferromagnetic or superconducting leads. In such a 
formalism, resonant transmission due to single particle 
interference, correlation effects arising from to strong 
electron-electron interactions, ferromagnetism and 
superconductivity proximity effect in the presence of 
ferromagnets and superconductors can be treated in a 
systematic way. We noticed that there exists a circuit 
theory for mesoscopic systems developed by Nazarov 
et al— based on the kinetic equations of quasiclassical 
Green functions^ which provides an alternative way to 
investigate the transport properties of hybrid structures 
with arbitrary connections. 20 However, such a formalism 
is not favorable to the systems of strong correlation, and 
apparently inapplicable to the cases where the single 
particle interference effect (for example, in resonant- 
tunneling structures) is prominent since the dependence 
on the relative coordinate of the quasiclassical Green 
functions is integrated out. 

Thanks to recent advances in nanofabrication and 
material growth technologies, several kinds of hybrid 
mesoscopic structures have been realized experimen- 
tally. These nanoscale structures include mesoscopic 
junctions such as normal-metal/superconductor— (N/S) 
and ferromagnet/superconductor 23 (F/S) contacts, 
superconductor/insulator/superconductor 24 (S/I/S) 
and superconductor/ferromagnet /superconductor— 
(S/F/S) junctions, and certain kinds of reso- 
nant structures such as superconductor/quantum- 
dot/superconductoi"2& (S — QD — S), normal- 
metal / superconducting quantum-dot /normal-metal— 
(N — SQD — N), normal-metal/ferromagnetic-quantum- 
dot/normal- metals (N — FQD — N) transistors. In a 



normal- metal/superconductor (N/S) junction, Andreev 
reflections^ dominates the transport process at low 
bias voltages, in which an electron in the normal metal 
slightly above the chemical potential of the supercon- 
ductor is reflected as a hole slightly below the chemical 
potential at the interface between the normal metal and 
superconductor with an electron pair moving into the su- 
perconductor, and vice versa. When two superconductor 
components are coupled together through an insulator 
or a normal metal, electron pairs can move coherently 
from one superconductor to the other, yielding a nonzero 
current even in the zero bias limit-the well known DC 
Josephson effect, and an oscillating current at finite 
bias-the AC Josephson effect— The S/N/S Josephson 
junction must be a mesoscopic system, with the length 
smaller than the phase coherent length of electrons in 
the normal region, to ensure electron (hole)'s coherent 
motion inside the normal part. Then a reflected electron 
(hole) can interfere constructively with itself, a process 
that produces a set of decoupled forward or backward 
'Andreev energy levels' carrying positive or negative 
Josephson current^ 2 ' 33 An impurity inside the normal 
region couples the Andreev energy levels, and thus 
modifies the quasi-particle energy spectrum and other 
quantities. In the presence of a ferromagnetic metal, 
a spin-polarized current may be induced due to the 
imbalance of the spin populations at the chemical 
potential.— Spin imbalance also introduces a net mag- 
netic moment-the magnetization of ferromagnets. When 
two ferromagnets participate in a transport experiment, 
the relative orientation of the magnetizations of these 
two ferromagnets will play an important role in the 
transport properties, and the spin-valve effect arises. 35 
Combining ferromagnets and superconductors, one may 
expect some new transport features, since there is no 
complete Andreev reflection at the F/S interface. The 
conductance of a F j S junction can be either smaller or 
larger than the N/S case, depending on the degree of 
spin polarization of the ferromagnet— When ferromag- 
nets, superconductors and confined (interacting) normal 
metals are integrated together, the interplay between 
ferromagnetism, superconductivity and electron-electron 
interaction is anticipated to lead to more interesting 
and more complicated transport properties. Despite the 
basic interest in the fundamental theory as mentioned 
above, hybrid mesoscopic systems also boast potential 
applications in future electronic devices which employ 
both the charge and spin degree of freedom of electrons. 

Starting from a microscopic Hamiltonian, we derive in 
this paper a general current formula within the Keldysh 
NEGF formalism for hybrid mesoscopic systems in which 
a central nanoscale interacting normal region is weakly 
connected to two leads, each of which is either a fer- 
romagnet or a superconductor, thus providing a uni- 
fied theory of electron transport in general hybrid struc- 
tures, which incorporates resonant tunneling, strong cor- 
relation, ferromagnetism and superconductivity proxim- 
ity effect. Such a formula can be also applied to the 
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nonequilibrium situation. Rather than from the original 
mean-field Stoner ferromagnet^ and BCS superconduc- 
tor Hamiltonian^ we calculate the current from their 
diagonalized forms after appropriate Bogoliubov trans- 
formations, with which the ferromagnetism and super- 
conducting proximity as well as the chemical potentials 
of the system are embodied in the tunneling parts of the 
system Hamiltonian. Such a procedure is found to be 
a crucial step in the analysis of our transport problem, 
and facilitates the applications of the general formula to 
the specific forms for given structures, which are Meir- 
Wingreen-type formulae.' 1 ' 2 Employing such a procedure 
it is very easy to check whether such a theory satis- 
fies the condition of gauge invariance, a requirement of 
all transport theories. Moreover, the energy-dependence 
and bias- voltage-dependence of the level-width functions 
and the distribution functions in the current formula 
are derived in a strict and natural way, while this has 
been done somewhat phenomenologically in the other 
formalisms Ai2iiSiiSiil This merit allows us to investigate 
the I-V characteristics of hybrid mesoscopic systems with 
a much more broad bias region. In addition we demon- 
strate that only their relative value for some physical 
quantities appears in the current formula after some uni- 
tary transformations. These quantities include the chem- 
ical potential, magnetization orientation of ferromagnet 
and the phase of the superconductor order parameter. 
Physically only their relative value can be measured in a 
transport experiment for these physical quantities, thus 
justifying the ad hoc assumption that one of them can 
be always set to zeiom^^^^ Such a formalism can be 
directly extended to the cases with more than two ex- 
ternal leads, which can be either ferromagnetic or su- 
perconducting. A shorter paper which summarizes the 
formulation has been reported elsewhere^ 

In order to illustrate the validity and versatility of our 
formulation, we apply the derived formulae to a non- 
interacting double-barrier resonant structure (DBTS) be- 
yond the wide-band approximation which is usually used 
in the Keldysh NEGF formalisms. 9 i 12 i 14 i 15 i 16 We neglect 
the interaction effects, since in a regime where these inter- 
actions are not important, we can then see more clearly 
how ferromagnetism and superconductivity influence the 
transport properties of a normal metal resonant struc- 
ture coupled to ferromagnetic and/or superconducting 
leads. As demonstrated in Section III, we derive the fi- 
nal current formula based on our formulation in a more 
systematic and economic way than others i^lL Some un- 
expected and novel transport features are found. When 
the level spacing of the central normal region is compara- 
ble to the bandwidth of the ferromagnetic lead(s), the I-V 
curves show resonant peaks plus resonant shoulders, re- 
flecting directly the profile of the density of states (DOS) 
of the Stoner ferromagnet. This observation provides an 
alternative way to measure the degree of spin polariza- 
tion of the system. The tunnel magnetoresistance (TMR) 
decreases non monotonically, as well as oscillates, as a 
function of the applied bias voltage between the ferro- 



magnetic leads. It is enhanced or suppressed within the 
resonant regions depending on the couplings to the two 
sides. We also find negative TMR at some bias voltages 
in the strong coupling limit. These features tell us that 
there is richer physics in the TMR of a resonant struc- 
ture. In the presence of ferromagnetic and superconuct- 
ing leads, a series of peaks emerges in the Andreev cur- 
rent whenever the resonant Andreev reflection condition 
at the N/S interface is satisfied as the gate voltage ap- 
plied to the central part varies. The number and height 
of these Andreev current peaks are strongly dependent 
on the bias voltage and the degree of spin polarization 
of the ferromagnet lead. Interesting step and peak struc- 
tures are observed in the I-V characteristics, which may 
be used to determine the DOS of both ferromagnetic and 
superconducting leads. Finally we investigate the DC 
Josephson current in S/I/N/I/S structures. It is shown 
that the DC Josephson current is slightly enhanced if the 
energy gaps of superconductors becomes asymmetric, in 
contrast to the S/N/S systems. The current-phase rela- 
tion is also weakly dependent on the asymmetry of the 
superconductor energy gaps. 

The rest of this paper is organized as follows: In Sec- 
tion II the full Hamiltonian of an interacting normal 
metal placed in between either ferromagnetic or super- 
conducting bulk leads is given. We first express the 
current in terms of the nonequilibrium Green functions 
in the Nambu-spinor space in general cases, and then 
present current expressions for the specific structures. 
The gauge invariance is proven to hold as well. Section 
III is devoted to the applications of the current formu- 
lae derived in Section II to non-interacting double-barrier 
structures, with a detailed analysis based on the analyt- 
ical results and numerical demonstrations. Concluding 
remarks are given in Section IV. An appendix is included 
to present the expressions of the self-energy matrices and 
level-width matrices due to the elastic couplings to the 
leads. 



II. FORMULATION OF THE PROBLEM 

We consider electron motion along the longitudinal di- 
rection i in a hybrid sandwich structure schematically 
shown in Fig. 1. The central part is assumed to be 
in the normal state, connected via tunnel barriers (in- 
sulators or point contacts, etc.) to two bulk materials 
acting as leads, each of which can be either a ferromag- 
net, or a conventional BCS superconductor. We adopt 
the Stoner model 3 ^ 3 ^ for the ferromagnet and the BCS 
Hamiltonian 31 for the superconductor. The Stoner model 
Hamiltonian is characterized by a mean-field exchange 
magnetization h, and can be written as 

H p = j dx*\x)(-^l-it-\L-ty{x), (1) 

where m* is the electron effective mass, a = {<? x ,a y ,a z ) 
is the Pauli spin operator, fi is the chemical potential 
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FIG. 1: A schematic diagram of a two-terminal hybrid meso- 
scopic structure. A mesoscopic normal region is attached to 
either ferromagnetic or superconducting leads with chemical 
potential fic and fin- In the former case the magnetization 
h makes an angle Of relative to the normal z axis. The cur- 
rent is assumed to flow along the longitudinal x axis from the 
higher-chemical-potential lead to the lower-chemical-potential 
one. 



and = (^j,^!) is the spinor field. In what follows 
we assume that the magnetization h makes an angle Of 
relative to the z axis, while we ignore the orientation 
with respect to the plane perpendicular to the transport 
direction being not relevant to the transport properties^ 
Within the mean-field approximation, the BCS Hamil- 
tonian takes the form 

U 

+ J dx[A{x)V\(x)y\(x) - A*(s)*i(a0*t(a0]-(2) 

In Eq. J2J) is the field operator of electrons with spin 
a =1,1, A(x) = t7(*f is the off-diagonal pair 

potential, with U a negative constant characterizing the 
electron-electron attraction. In general the pair poten- 
tial A needs to be determined self-consistently, and in 
this work will be assumed to be position and energy in- 
dependent for simplicity. 



Since we are concerning about transport properties 
of electrons rather than their motion in real space it 
is more convenient to deal with the lead Hamiltonian 
in fc-space. Expanding the electron field operator in 
terms of the eigenfunctions of the momentum opera- 
tor in the longitudinal(tunncling) direction as "to- (a;) — 

fkae lkx (ska-e tkx ) we cast the ferromagnet and super- 
conductor Hamiltonian and J3J) into the following 
forms in fc-space (subscript 7 = C and 1Z are added to 
denote which side of the structure the ferromagnet or 
superconductor is located at) 



n (F) 



X^[ £ 7fc - sgn(a)h-y cos0 7/ - /x^/j^/yfc, 

feer 

ka 

E^ ~ M 7 ) S 7 fc<T S 7fe^ + 

fe<7 

E { A ~r s lkV s l-H + A 7 S 7 -fc|S 7fet , (4) 



where e 7 & = h 2 k 2 /2m*, a stands for the opposite to a. 
Here the order parameter A 7 is characterized by its mag- 



nitude and phase: |A-, 



and, as we will show, |A„ 



opens an energy gap in the excitation spectrum of the 
superconductor. / 7 fca-(/ 7 fc CT ) and s 7fecr (s^ fco .) are the elec- 
tron destruction (creation) operators of spin a in state 
k in the ferromagnet and superconductor, respectively. 
Henceforth, the notations a =| , J, and a = ± are used in- 
terchangeably. In what follows physical quantities such 
as the particle operator ip 7 particle energy e, and the 
chemical potential /i of the different parts are labelled by 
subscript 7 = £,1Z;C, whenever convenient. In addition, 
ferromagnetic or superconducting characteristics of the 
leads is stressed by adding a subscript or a superscript / 
or s to some quantities. 

The Hamiltonian of the central region Tic m momen- 
tum space can be modeled by 



1 la 

(5) 



where ^] no .(?/> crlo -) creates (destroys) an electron of spin 
a in state n, and TLi n t represents the interaction terms in 
the central region. It may include the electron-electron 
Coulomb interaction 



H 



d — el 



E 



U: 



na, ma' 



IpLatPcnv^l^lpcna' , (6) 



n, m, a, u 
na ma' 



or the electron-phonon interaction 



H 



el—ph _ ^""^ 
9 



fiUqCltq + E U na,q^l na ^cna{C,l +(-q)- 



n,a;q 



(7) 
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Here the first term is the free-phonon Hamiltonian, while 
the second represents the electron-phonon interaction, 
with interaction matrix element U na>q . Q(C«) is the 
phonon creation (destruction) operator in mode q. 

The couplings between the leads and the central re- 
gion can be modeled by tunneling Hamiltonian, no mat- 
ter how the leads are coupled to the central region, pro- 
vided the couplings are not strong enough. Certainly the 
coupling strengths depend on the detailed configuration 
of the setup and should be determined in a self-consistent 
manner. However, for simplicity they are assumed known 
and can be written as 



kn\a 



(9) 



kn:a 



To see tunneling processes more clearly, and, more im- 
portantly, to facilitate the analysis of gauge invariancc 
and the simplification of the general current formula (|3 1|) 
to the forms of specific systems, we first diagonalize the 
Hamiltonian of the leads by Bogoliubov transformations. 
For the ferromagnetic lead one has 



f lka = cos(0 7 //2)?/> 7 /fc<x - sgn(a) sin(0 7 //2)V> 



7/k<7 



(10) 

and for the superconducting lead, 
e'^^Sjka = cos9 lsk ^ lska + sgn(a) s,inQ 18k Vi^ _■ 

*ysk(7 
(11) 

In Eqs. p()lll|l . V*(V) is the pair creation (destruction) 
operator guaranteeing particle conservation, which trans- 



forms a given N-particle state into an (N+2)-particle ((N- 
2)-particle) state, i.e., /V\N) = \N + 2)/\N -2), and 



9-ysk = arctan 



'e^. + ^e 2 lk + \A 1 \^ 



1/2 



(12) 



Substituting the Bogoliubov transformations l|lUfl and 
(111(1 into the lead Hamiltonian © and Q, we get the 
following diagonalized forms for the ferromagnetic and 
superconducting leads, respectively, 

ka 

= ^Z e lfk<y^l fka ^fka, (13) 



w f = Y,N^ k + \^~^ lsk ^ lska + 



ka 



Constant 

= E £ i sk ^\skai ) iska + Constant (14) 



ka 



Now the particle number operator commutes with the 
corresponding lead Hamiltonian. After the Bogoliubov 
transformations, the diagonalized lead Hamiltonian de- 
scribe the excitation (quasi particle) properties. The 
minimum energy of the excitations in a superconductor is 
|A 7 |, implying an energy gap in the excitation spectrum. 

With the Bogoliubov transformations (|l()|l and (|11|1 we 
turn the tunneling Hamiltonian into 



n 



7 (F) 



E V kl,„ [cos(07//2)V> 7/fe(T - sgn(a) sin(fl 7/ /2)^ _] + 

7 /fco- 



V kn*a^lna [ cos(6> 7/ /2)V> 7 /fc CT - sgn(a) sin(0 7 //2)^ 



7/fecr- 



Hp S) = E ^[cos^^ + ^n^sin^^ ^pt] gW*^ 



KT>L,e-W 2 [cos8 lsk ^ ska + sgn(a) sm(9 Jsk )P^ 



^jska 



r 



(15) 



(16) 



The associated physical processes are more obvious and 
clear in the semiconductor modeli^^i an electron of spin 
a in the central regime can tunnel into cither the spin a 

band or a band of the ferromagnetic lead, or tunnel into 
a spin a state or condensate into an electron pair with a 
hole state of opposite spin being created, and vice versa. 
In superconductors, correlation between two creation 



or annihilation quasi-particle operators with opposite 
spins are very important, relating to the Andreev re- 
flection in transport processes. When ferromagnets are 
introduced, the correlation between a creation and an 
annihilation quasi-particle operator with opposite spins 
needs to be considered. To incorporate these two kinds 
of correlations in a unified way and consider the ferro- 
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magnet and superconductor on the same footing, we here 
introduce a generalized Nambu-spin representation span- 
ning a 4-dimensional spin-orientated particle-hole space 

= Ol T ipxl 4>li Vm) 1 - Within the Keldysh 
NEGF formalism, Green functions are defined as 

G x , y (t 1 ,t 2 )=i{Tc[& x (t 1 )®*l(t2)}), (17) 

where Tq is the time-ordering operator along the 
closed time pathiS The usual retarded/advanced and 
lesser/greater Green functions then take the form 

= T«(±tiT*a)X)<[*o*(*i)® 

*J J -(*2)+*J J -(*a)®*ai(«l)]), 

G^(h,h) = EG</>( tl ,t 2 ) 

•J 

= ±^(*«(*2)/*ai(*l)® 



in which 



^(tx)/*^.^)), 

where a, /? = 7/, 7s, c and i, j = fc, n. 

In this 4-dimcnsional Nambu-spinor space the total 
Hamiltonian can be rewritten in the following compact 
form 



TC = Tic + Tic + 7~tn + Ti.T + 'Ht, 



(18) 



where 



H 



He = ^fLE cn * cn + 

n 

H mt ({*in,*cn}), (19) 

€ /S) = £*; //sfc E 7/M * 7/M , (20) 

fc 

= EK//X /s c)*™+^4 (21) 



7 

7(F/S) 

r 



fcn 













































— e^T 



In writing down Eq. (|19I20I21|1 . we have introduced the 
energy matrices 



E, 



and the tunneling matrices 

V#(*) = R / (^f)V^P( M7C 0, (23) 

V£(t) = R'(0 7 dk)V&P(^t+^), (24) 
M7C = M7 _ 




V 



7/A 
kn 



R/(a 



R%x) = 



P(z) 






V 

/ cos a: 


- since 

V 

/ cos a: 
V*sinx 







-v^<" 





Tills 



7/A 



V 

/ e fa / fi 

e -«/fi 







sin a; 

cos x — sin x 

cos x 

sin x cos x 

— Psinx 
cos a; 
cos a; 

-"P*sinx 







t/7//s* / 
V kn / 






■psina; 
cos a; 





V 






ix/h 











are the coupling, rotation and phase matrices, re- 
spectively. Note that we have performed a gauge 
transformation's to get the above Hamiltonian. Chem- 
ical potential is now incorporated in the phase opera- 
tor P(/i 7 ct) in the tunneling matrices Eas. i|23l24|) . which 
along with the rotation operators is very useful to demon- 
strate gauge invariance for our system as shown below. 

The current flowing from lead 7 — C, 1Z to the cen- 
tral region can be defined as the rate of change of 
the electron number iV 7 = J2 ka flkaf"r^(s\ k(T s 7krT ) = 

Y.ka^lfka^ifkMlska^iskv) in the lead. Within the 
Keldysh NEGF formalism, the current is expressed as 

l y (t) = -e < 7V 7 >= — < [N y ,Hr] > 

([G<, 7//sfe (M)V^ /s (i)- 

yii /s Ht)G< f/sk>cn (t,tt 

{(vL /M (*)G< M>on (M))J. 

(25) 



i— 1,3 nk 



nk 



Since the Hamiltonian of lead 7 is of the form 
^If/ska^if/ska, the equations of motion for G< f/skcn 
along with the Langreth analytic continuation^ yield the 
following Dyson equations 



Q< 

7/ j sk.cn \ 



G 



< 



(*!,*') +g 7 C /M , 7/M (*,*l) 



V fcm S (*l) G cm,cn(*l; 



(26) 



cn.,7// sfc 



cn.cm 



(*.*i)V& M (ti) 



7 



g: 



V^ M (ii)< //sfc , 7/M ^i^')j,(27) 

in which the unperturbed retarded/advanced Green func- 
tion Sjf/ a 1 f/ s °f l ea< i 7 can De readily obtained from the 
Hamiltonian l|13l 14|l as diagonal matrices 



S 7 // S fc, 7 // S fe (M') 











o \ 



















9'jf/sk 





V 








T+ 
9 7 /M / 



(28) 



(29) 



and the lesser (greater) Green functions are re- 
lated to the retarded (advanced) Green func- 
tions by gj/£ fc , 7/M (i,t / ) = [f 7 ( e 7/M) - ± 



2-*-] [s 7 //sfe, 7 //sfc(*>*') s 7 // s fe i7 // s fe Cm')] ■ 

distribution matrix f 7 (e 7 y/ s fc) reads 

f 7 (e 7 //sfc) = 
/ /(£ 7 /Mt) 

o /(- 






The Fermi 



V 



o 

r 7/Ml) 









/( £ 7/Mi) 

/(-£ 7 // s fet). 



where f(x) = (1 + e x / kBT ) 1 and we have used the rela- 
tion f(-x) = 1 - f(x). 

Substituting Eq. 126() into 1|25|1. we obtain 



2 7 (t) 



i=1.3 ,i 

ze x - 
T 2^ 



*iRe|^ 



1 cm.cn 



7 //s;nro(*> *l) 

(*i,t) + £< //s . nm (i,fi) 

Gcm,cn(*l, *))..}, (30) 



where 



s^c^) = E v ^ /St (^) 



r/a,</> , \\fyf/ s (+ \ 

^f/sk, 1 f/sk\ t ^ t ^> y km y t2 > 

is the self-energy matrix (see Appendix A) arising from 
electron tunneling between the central region and lead 7. 

For steady transport, no charge piles up in the central 
normal region. One then has 2x(i) = —In(t)m& After 
symmetrizing the current formula (|30|l . we finally get by 
setting l(t) = [l c (t) - I K (t)]/2 



l(t) 



i=l,3 



,(Mi 



^ //4 ,(t,*i)]G< c (t ls t)+[S< //a (t,t 1 )- 
E< //a (Mi)]G£ (ti,t)) }, (31) 



where the trace is over the level indices of the central 
region. Eq. (|31(1 along with the self-energy matrices 
given in Appendix A is the central result of this work. 
The current is expressed in terms of the local properties 
(G r / a ) and the occupation (G < / > ) of the central inter- 
acting region and the equilibrium properties (S < ^ > ) of 
the leads. It is emphasized that the current is usually 
independent of time except the presence of two supercon- 
ductor leads with nonzero bias voltage. Formula (|31|l can 
be employed to investigate both equilibrium and nonequi- 
librium electronic transport in various kinds of hybrid 
mesoscopic systems, including F/I/N/I/F, F/I/N/I / S, 
F/I/N/I/N, S/I/N/I/S, S/I/N/I/N, and F-QD-F, 
F-QD-S, F-QD-N, S-QD-S, S-QD- N struc- 
tures as well, in which arbitrary interactions are allowed 
in the central part of the structure. 

It is not difficult to check that Eq. I|31[l is gauge in- 
variant, i.e., the current T(t) remains unchanged under 
a global energy shift in the whole region. This can be 
achieved through a gauge transformation for the Hamil- 
tonian of the system 



f(e i) = exp j^eotC^VLT^c 



E ^If/sk^-yf/ska)}, 



j—C,TZ;ka 

where eo is just the energy shift. Such a gauge transfor- 
mation gives rise, in turn, to the phase transformations 
of all the terms in the right-hand side of Eq. (|31[) 

s;/< s (Mi) - P(eo*)E;/< (MOPtfatx), 

G a c /<(t u t) - pt^OG^^PM)- 

The above procedures, equivalent to applying a 
phase transformation to the current operator: 
T 3 ' (eot)T(t)T?(eot), ensure that the current remains 
the same under such a transformation. Therefore, the 
current formula l|31|) is gauge invariant. 

Now we check whether the current becomes zero if we 
take the zero bias limit /i£ = fi-jz — fJ-o- We first perform 
a phase operation P(//oe^ + ^k/2) corresponding to the 
gauge transformation f(/j,cot), /ioc = — to Ea. (|31|l . 
obtaining 

iw--jE/s"^{(5[^"w-^"w] 



G<M - [f'f'MWt) - r*'"(t)f K ( e ) 



where 



7e 

G a cjs))\, 



(32) 



t cf/ S{£) = pT( ^ )r «/. (e) p ( p > 



g^/<(M)pW + ^), 
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with if s = (p£—(pn- From the fluctuation-dissipation the- 



orem G c ' c 



[f e9 (e)-il±il](G« -G-) (f £ =fv 



1i 



f eg ), one can readily verify that the current is zero except 
in the presence of two superconductor leads with different 
superconducting order parameter phases. In this case, 
there still exists a DC Josephson current in the zero bias 
limit due to the coherent tunneling of quasi-particle pairs. 
This can be seen more clearly in the expressions of the 
current in the specific systems IJ33I3 



Up to now we have obtained the expression for the 
current in a general case in which each of the two leads 
can be either a ferromagnet or a superconductor. Next 
we apply this general result (|31|1 to the specific structures 
we are interested in. We first consider the case in which 
two leads are ferromagnetic. Inserting the expressions 
of the self-energy matrices X 7 /(ii, t^} (Appendix A) into 
Eq.lEU, we get the current in a F/I/N/I/F or F- QD- 
F structure after a rotation transformation and a phase 
transformation 



i=i,3 



G< c (e)+[f c f(sTeV)i c (eTeV)- 
T^(e)i K (e)][Gl c (e)-Gl c (e)]) }, (33) 

where f c ? = R/t(^-)r £ -f R/(^), 8 f = Cf - 6 nf , and 



£v(e=Fc) 



(f{e-c) 

f(e + c) 

f{e-c) 

V f(e + c), 

The expression of tunneling current l|33|) resembles 
formally the current formula derived by Meir and 
Wingreeni^ for a confined region coupled to two normal 
electrodes. The difference lies in that the coupling ma- 
trices and Green functions in (|33l) arc spanned in the 
Nambu-spinor space, which reflects the dependence of 
the current on the spin polarization of the ferromagnetic 
leads and the relative orientation of the magnetic mo- 
ments. When we set to zero the magnetic moments of 
the two leads, Eq. l|33l reduces to Equation (5) in the 
paper of Meir and WingreenjiS since in this case the fer- 
romagnetic leads become normal metals. As seen from 
Eq. (|33|l . current is only dependent on the relative ori- 
entation of the magnetizations of two leads, although 
there is an apparent 9-jif dependence in the expression 
for G r > a /<. Nevertheless, this dependence of the Green 
functions G r > a / < on the orientation of the ferromagnet 
magnetization comes from the self-energy matrices T, ' 
and Y, n f after the rotation operation R^(^ i ), hence 



they only depend on the relative orientation as can be 
seen more clearly in the non-interacting model. 

If one lead (£) is ferromagnetic and the other (72.) is 
superconducting it is expected that Andreev reflection 
process, dependent on the spin polarization of the ferro- 
magnet, will dominate the current at low bias voltages. 
Applying a phase and a rotation transformations to Eq. 
(|31|1 . simple integration gives 

G< c (e)+[T c f(eTeV)ic(eTeV)- 
Tf (e){ K (e)} [G; C ( £ )-G° C ( £ )])J, (34) 

in which the full Green functions are 

Viz, 



= J d{t-t l y^ t - t 'v h v( mc t + ^-) 



R/(^)G-/<(M')R /t ( : f)P t (M^t'+^). 

In Eq. H34fl the current does not depend on the orienta- 
tion of the magnetization of the ferromagnetic lead and 
the phase of the order parameter of the superconductor 
lead. This can be clearly demonstrated by expanding the 
full Green functions of the central part perturbatively, as 
we will show below in the non-interacting case. In ad- 
dition, one can divide the current into several parts im- 
plying the contributions from different physical processes 
such as normal particle tunneling and Andreev reflection, 
after expanding the right hand side of Eq. i|34[l . We will 
show it later in the non-interaction case. 

When two leads are superconducting, the situation be- 
comes much more complicated. As did in the previous 
examples, we derive the following current formula for 
S/I/N/I/S or S - QD - S systems 



'Fsns (^) 



1=1,3 

e \ ^ 



Tf s (s)]G< c (e;t)-[tf s (e T eV;t) 
i c (sTeV)-rf s (e)i n (e)]G a c , c &t)) }, (35) 



where 



T c e J p (sTeV;t) = ^(eVt + ^)T%(e T eV) 



P(eVt + 



G r c <;J<(e;t) = fd(t-t')e i ^- t 'y n P( m et+^) 

g^/<(m')pW*' + ^), 

with ip s — ipc — fiz- One may wonder why we use 
the notation G^^ (e;t) with the additional variable t 
other than Gc,c^ < (e). The reason is that the full Green 
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functions Gc,c should be calculated in the presence of 
tunneling between the central part and the two sides, 
as well as the interactions in the central region. In the 
present case, the ^-dependence can not be avoided in the 
self-energy matrices, while it can be removed by a uni- 
tary phase operation when only one superconductor is 
involved. The current through a confined interacting re- 
gion connected to two superconductor leads is generally 
time dependent, as in the case of biased weak Josephson 
linksi^ii However, in the limiting case of zero bias, the 
current is a time-independent nonzero quantity, as can be 
seen from Eq. (|35|l . In other theoretical treatments 
a double-time Fourier transformation is usually taken as 
X(t,t') = ±J2 n I duje- luJt e l ^ +nuJ °/ 2 '> t ' 'X{u,u + nu} /2), 
where ujq = 2eV/h, and the current yields a har- 
monic expansion of the fundamental frequency T(t) — 
En^e™" '. In fact, the Green functions G^ 7< in J33 
can be expanded in powers of the fundamental frequency 
co, i.e., G^M) = E m G^ /< (e, £ +mc /2) e — ot/2 ; 

which with the expression for the Green function G r J a J < 
below Eq. 1|35|) is exactly of the form of the double-energy 
transformation. 42 However, we show here that one can 
obtain in principle the time dependence of the current, 
as long as one can derive the full Green functions of the 
central part, which need further investigation. 

So far we have presented a general formulation to cal- 
culate the current through a confined normal region con- 
nected to two leads being either ferromagnetic or super- 
conducting, and the Meir-Wingreen-type formulae in the 
specific cases. Although the formalism can not applied 
to the strong coupling of the central normal part to the 
outer world as the circuit theory of the hybrid meso- 
scopic transport^ it permits us to investigate the effects 
of the single particle interference and strong electron- 
electron interaction on the transport properties of hy- 
brid mesoscopic systems, which is ignored in the circuit 
theory. Compared to the other formalisms based on the 
similar Keldysh NEGF technique t 9 i 12 i 15 i 16 i 17 i 67 our for- 
malism is more systematic and more general. In the 
present formalism, one do not need to make additional 
ad hoc assumptions as mentioned in the introduction. 
One can also judge what quantities can be measured in 
experiment in a much more complicated structure, by 
observing simply the energy-independent arguments in 
the exponential functions or the triangle functions in the 
unitary matrices of the tunneling parts of the Hamilto- 
nian after the Bogoliubov transformations. As we will 
demonstrate, we obtain the final current formula after 
simple matrix algebra rather than difficult mathemati- 
cal techniques pi^iLiLgi This kind of mathematical sim- 
plicity makes our formalism more appearing than others. 



More importantly, an explicit energy- and bias-voltage 
dependence of the level-width functions and distribution 
functions allows us to investigate the I-V characteristics 
in a much more wider range of bias voltage. 

III. APPLICATIONS TO THE 
NONINTERACTING MODEL IN THE CENTRAL 
REGION 

In this section we use the formulae developed above 
to study transport properties of various kinds of hybrid 
mesoscopic systems in which, for simplicity and conve- 
nience of comparison with other theories, the interac- 
tion effects in the central confined region are not consid- 
ered. The absence of the interactions permits an anal- 
ysis of the genuine physical influence of ferromagnetism 
and superconductor proximity on the transport proper- 
ties in hybrid structures. One of the best candidates for 
a non-interaction confined region is double-barrier res- 
onant structures (DBTSs) with quantized discrete en- 
ergy levels^ Therefore we adopt the double-barrier model 
with the emitter and collector replaced by either a ferro- 
magnet or a superconductor. Throughout the following 
calculations we use the following approximations: (i) the 
level shift is omitted, (ii) the coupling coefficients are 
real constants, independent of spin and energy such that 
V^lj = V~<il s = yifls*. However, we will abandon the 
ugua y9 i i2 i i4 i i5 i i6 w ide-band approximation of the level- 
width functions which is reasonable in the low-voltage 
transport, since in this paper we also deal with high bias 
voltage situation. As we will show, this permits us to 
investigate the current within a much wider region of 
bias voltage and find interesting transport features of 
the same resonant structures which can not be found 
in the other formalisms based on the Keldysh NEGF 
technique P iiiSi 1 i 

In the absence of interactions within the intermediate 
normal metal, the full retarded/advanced Green function 
can be solved from Dyson equation 

G% a (t,f) = G Q ;J a (tX) + J dt x J dt 2 G°;;/«(Mi) 

V r / a (h,t 2 )G r c / r a (t 2 ,t') 
= G^ a (t,t') + J dt x J dt 2 G r J«{tM) 

^ a (h,t 2 )G^ a (t 2 ,l/), (36) 

in which G c , c is the decoupled Green function, which 
becomes when the central region is isolated from the out- 
side world 
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-tl)l% Q 

e i(e„i-fic)(,ti-t2)/fi 









-Mc)(*i 




-t 2 )/ft 







,T-A«c)(ti 



(37) 



The lesser/greater Green function of the central region 
is calculated via Keldysh equation 

G</ > (f,t / ) = (dtx J dt 2 J dt 3 J dh 

[l + G r CtC {t,t x )^ r {tiM)] 
G°< /> (i 2 ,t 3 ) 

[l + £ (t3,f 4 )G« c (i 4) t')] + 
J dt x J dt 2 G r CyC (t,ti) 

£</>(*!, ta)G» c (t 3l O- (38) 

Once the full Green functions in the central region are 
known, we then have the complete knowledge to inves- 
tigate tunneling processes in the specific structures. In 
the following we study electron tunneling in three typical 
hybrid structures: (a) F/I/N/I/F magnetic DBTSs, (b) 
F/I/N/I/S DBTSs, and (c) S/I/N/I/S DBTSs. 

A. F/I/N/I/F structures 

When two ferromagnets are separated by a thin non- 
magnetic barrier, two kinds of physical effects arise. The 
first is the spin valve effect showing a (1 + ecos6*) 
dependence of the tunnel conductance on the relative ori- 
entation 9 between the involved two magnetizations. The 
other is tunnel magnetoresistance (TMR)j^ 5 . defined by 
AR/R = (R a — R p )/R a , where R p and R a are the re- 
sistances when two magnetizations are parallel and an- 
tiparallel, respectively. The spin valve and TMR are 
due to the spin polarization induced by an exchang- 
ing coupling between electron spins and the internal 
magnetizationfi&iS and the relative orientation of the 
magnetizations can be adjusted by applying a magnetic 
field. A TMR up to 11.8% at room temperatures, and 
24% at A.2K, was reported in CoFe/Al 2 3 /Co planar 
magnetic junctionsi-^ It is observed that the TMR de- 
creases with increasing bias voltage^ 

Recently, double-barrier magnetic resonant struc- 
tures has attracted much experimental 44 and 
theoretical 16 i 45 i 46 i 47 i 48 i 49 attention. The theoretical 
results show that the TMR of resonant magnetic 
structures is enhanced compared to the single mag- 
netic junctions due to resonant tunnelingi^^ii 7 ^ 8 ^ 9 . 
However, there appears a controversial issue related 
to the TMR enhancement, on whether it is for peak^ 5 - 




FIG. 2: A schematic potential profile for a biased magnetic 
double-barrier structure with two ferromagnetic electrodes 
characterized by their respective magnetizations he and hn- 
The hatched regions denote the states occupied by electrons. 



or valley currenti* 6 - In addition, Sheng et all* 7 - found 
both positive and negative TMR in F/I/F/I/F double 
junctions. To clarify these issues, we re-examine this 
problem using our formula (|33JI in terms of the two-band 
free-electron spin-polarization models for the ferromag- 
netic leads. A typical double-barrier magnetic structure 
is schematically plotted in Fig. 2. 

We take two steps to calculate the full re- 
tarded/advanced Green function G r J.c{t, t') of the central 
region. First we decouple the system from the left ferro- 
magnet, denoting the corresponding retarded/advanced 

Green function by {7c,c°(£, i')j then we couple the cen- 
tral region to the left ferromagnet and calculate the full 
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Green function G from Q. From the Dyson equation OBI Substituting the self-energy matrices E~ ? (Appendix A) 



we have 

g^(m') - £/°(t,0 + / *i / dt 2 g r c / c a (t, tl ) 

TTjffa^GtffaJ), (39) 

%i a } {txM)Q:i c a {h,t'), (40) 

where (X = G r < a /<, £ r < a /<, g r < Q /< and 57> a /<) 

X(M') = P( MKC t)R/(^)X(t,t')R /t (^)P t (^ci')- 



into the Dyson equation l|4()|l. one has 



de 



(41) 



where the retarded/advanced Green function for the iso- 
lated central region is 



/(E — 



§#00 



(E 



o 








(E e -e' ±i0+ ) 





(E 









(42) 



with d na = e n(T — /ic + fin. The full retarded/advanced 
Green function is obtained in the same way 



G^(M') 



de 
2V 



-ie(t-t')/h 



df_ fe (i_t')/fi 
2^ 



~f^( £T eV) 



r^( £ )±^f^( £T ey)) 



/ ge-^-^G#( £ ). 



(43) 



The lesser Green function is associated with the re- 
tarded and advanced Green functions via the Keldysh 
equation i|38|) where g < (i,f') = 

G< c (t,f) = y dtx y dtaG^eC*,*!) [S<f(t!, t 2 ) 

+S< / / > (t 1 ,t 2 )]G^ c (t 2 ,t') 

= y g e ^ (t - t,)/?l G^( e )[f^( e Te^) 

ic{eTeV)-T^f{e)in{e)]Gl c {e) 
= J ^e- fe (*-*')/ R G< c ( £ ). (44) 

Substituting the advanced and lesser Green functions 
and (glj into Eq. or (J3TJ , one gets the following 



Landauer-Buttiker-type^ formula for the current through 
a non-interacting F/I/N/I/F magnetic structure, 

Zfnf(0f) = ^J de[f c (e - eV) - Me)]T /nf (e,9 f ), 

(45) 

where the transmission coefficient Tfnf is cast into the 
following compact form 

Tfnf(e,9 f ) 

i=l,3 

= \ [(cos 2 ^ff rf + sin 2 |f f rf )|^, c;11 | 2 



(2rf f rf / + cos^ ^r^rf 



|G; c;13 | 2 -(2r^r^ + cos 2 ^ff^ + 
sin 2 £ff r^)|G; c;31 | 2 + (cos 2 ^ff rf + 



sin 2 ^f^r^)|G; c;3 3l 2 



(46) 



Notice that we have dropped the arguments e — eV in f jr* 
and e in for brevity. The full retarded Green func- 
tion in Eq. (|46() is determined by the matrix inversion 

Gc, c = [i^c 1 + h^ Cf + i rTC/ ] _1 - and we have 

n "J- 
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sin 2 ^ff' 



-(cos 2 yff 7 -r -u,. — , 



-|sin^(rf -rf) 



Gr 
c,c;13 



= G 



c,c:31 



G 



c.c;33 



1 



e-el 



l(cos 2 ^ff +sin 2 ^ff + 1?' 



ff\ t ) — Lt c,c;11 Lj c : c;33 U c,c;13^ r c,c;31 • 

To obtain the last equality of Eq. (|46|1 . we have used 



[9 c ,c;ll ± o 1 llj^c, 



c;13 



fo r/a ,", l ±-r{,)G r/a , 1 = 



_|_ * p/ />/a 
^ 2 13 c, c;33 

_i_ 1 p/ /V/a 
31 2 31 "cjC;!!) 



where = +r TC/ . 

One sees from Eqs. lj4"5)l and (|4"rj)l that the current has 
a generic dependence on the relative orientation Of be- 
tween the two magnetizations. By observing the current 
expression and scrutinizing the structure of the Green 
functions, it is not difficult to find that the tunneling 
current through the magnetic structure is generally max- 
imized at Of = (parallel magnetization) and minimized 
at Of = 7r (antiparallel magnetization), a typical spin 
valve effect also in magnetic resonant tunneling devices 
(data not shown here). The ferromagnetism is reflected 
in the Of- and T^-dependence of the full Green func- 
tions of the central part G CiC , as well the transmission 
function T . When at least one lead is nonmagnetic, the 
^/-dependence can be removed with the help of a rota- 
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FIG. 3: I-V characteristics of a magnetic double-barrier struc- 
ture for (a) different magnetizations he and h-jz and (b) dif- 
ferent couplings r £ ^(0) and (0) at temperature ksT = 
0.1W. The level spacing is 1.5W, larger than the bandwidth 
W. 



tion operation R/. If we set Tf f = Tf f and Tf } = rf f , 
the current formula 1451) will recover the usual Landauer- 
Biittiker formula and the transmission is finally simplified 
to the Breit-Wigner type in the single level case. We no- 
tice that the current formula l|45|l is formally similar to 
the results of Wang et and Zhu et al.ji£ however, 
the discrepancy is nontrivial. The current formula 145(1 
allows us to calculate the I-V curves in a much wider re- 
gion of bias voltage yielding rich physics, while according 
to the theoretical treatments of Wang et alii& and Zhu et 
al. jii the current formula is restricted to the low bias volt- 
age case where the level-width functions can be viewed 
as energy-independent constants, and thus may result in 
even wrong consequences in the large bias voltage limit. 

Of particular interest is the current-voltage character- 
istics of double-barrier structures. In subsequent calcu- 
lations, we approximate the density of states of the fer- 
romagnetic leads by that of the two-band free-electron 
spin-polarization model 38,46 and take into consideration 
the finite width of these two bands. In this model 
the dimensionless DOS of the spin bands is pj^(s) oc 
y/(e + ah-y + W) /W, where W is the bandwidth mea- width is about (W — h 7 )/2h y . Moreover, the valley cur- 
sured from the band bottom to the Fermi level. For rent in a normal resonant tunneling structure is lifted 
some ferromagnetic metals, this is quite an appropriate when the leads become ferromagnetic. These surpris- 
ing results, unexpected within the wide-band approxi- 



in the central region, e n f = e„i = e n . To model the bias 
voltage drop inside the well we take /ic = fJ-n — eV/2, 
since the bias potential is assumed distributing uniformly 
across the double-barrier structure. Without any loss 
of generality we consider two quasi-stationary levels in 
the well, of which the energy of the lowest one is 0.25IT 
when V — and the level spacing is chosen as 1.514 7 , 
larger than the bandwidth of the ferromagnets. This as- 
sumption is, in practice, quite reasonable for the narrow- 
band ferromagnetic metals and quantum wells with very 
large level spacing or small quantum dots with very large 
charging energy. The I-V curves are shown in Fig. 3. 
The dotted line in Fig. 3(a), the well-known I-V char- 
acteristics of usual double-barrier structures, is given for 
comparison with the ferromagnetic case. In the pres- 
ence of ferromagnets, the structure of resonant shoul- 
ders neighboring to resonant peaks is celebrated in the 
I-V plots. The ratio of the peak width to the shoulder 



approximation 



38,50,51 



In the absence of a magnetic field 
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mation, can be understood from the potential profile of 
this kind of magnetic double junction structure, shown 
in Fig. 2. ft is well known that the current through 
a usual double-barrier structure is resonantly enhanced 
when one of the well levels falls into the region of a 
Fermi sea, i.e., eV < e n + eV/2 < W + eV& In the 
ferromagnetic situation, the Fermi sea is distorted and 
comprises two distinct parts: one with both spin-up and 
spin-down bands, and the other with only a spin-up or 
spin-down band, as displayed in Fig. 2. We thus have 
two types of Fermi sea for ferromagnets, one is repre- 
sented by h-y + eV < e n + eV/2 <W + eV and the other 
is — hry + eV < e n + eV/2 < /i 7 + eV. It is obvious that 
the resonant current through one of the well levels be- 
ing within the sea of the first type is larger than that 
of the second type, which is clearly reflected in the I-V 
characteristics in Fig. 3(a). 

Next we investigate the influence of coupling asymme- 
try on the tunneling current. The results are presented 
in Fig. 3(6). The magnitude of resonant current is sig- 
nificantly enhanced when one of the couplings becomes 
10 times as large. The coupling asymmetry induces a 
more significant enhancement of the tunneling current 
if it is the coupling to the higher-voltage lead (emitter) 
that is stronger, consistent with the tight-binding numer- 
ical result in the usual DBTSS The ferromagnetic I-V 
characteristics (peaks plus shoulders) in the reverse-bias 
case is blurred, also due to the same coupling asymme- 
try effect. These features can be understood in a similar 
way. The resonant current is roughly proportional to the 
ratio T c f{0)p c fT n f{0)p n f /{T c f{0)p c f + T n f(0)p n f) 2 , 
which becomes larger when one of the couplings T c ^(0) or 
r Kj (0) is enhanced. However, the magnitude of this en- 
hancement also depends on the DOS of the ferromagnet 
lead p^ f . If one strengthens the coupling to the lower- 
voltage lead (collector), the tunneling current is slightly 
enhanced since the DOS of the collector is comparatively 
large. 

Following Sheng et al.j 4 ^ we define the tunnel 
magnetoresistance (TMR) as AR/R = [lf n f(n) — 
2 "/™/(0)]/max(X / „ / (7r),X / „ / (0)). In Fig. 4 we give the 
TMR as a function of the bias voltage for some typical 
couplings, and the I — V curves in the cases of paral- 
lel and antiparallcl alignments of magnetizations for the 
convenience of comparison and analysis. In contrary to 
the monotonous decay with the bias voltage^ the TMR 
in magnetic DBTSs displays complex dependence on the 
bias voltage no matter what the values of couplings are, 
which arises from the resonant tunneling of electrons. 
This feature reveals that there is richer physics in the 
TMR of magnetic resonant structures. We notice also 
that the bias voltage dependence of the TMR can be com- 
paratively simple if the collector ferromagnet is of low de- 
gree of spin polarization h 1 /W , as shown in the negative 
bias domains of Fig. 4. This phenomenon can be ascribed 
to the weak perturbation of the spin-up and spin-down 
DOS of the ferromagnet with small h 1 jW by an external 
magnetic field. In addition, the peculiar behavior of the 



TMR also depends on the strengths and symmetry of the 
elastic couplings. In the strong coupling case (Fig. 4(d)), 
the TMR shows a resonant behavior similar to that in the 
tunneling current^^ and can even be negative at some 
bias voltages4i In other cases (Fig. 4(et)-(c)), the TMR 
drops within the resonant peak region and then develops 
peaks at the boundaries between the current peaks and 
shoulders, similar to the result of a non-interacting quan- 
tum dot coupled to two magnetic leads4& The different 
oscillatory behaviors of the TMR with the bias voltage for 
different couplings imply that the analysis on the TMR 
in the resonant structure 4 ^ may not stand. It is worth 
noticing that the TMR will eventually decay to zero in 
the large bias voltage limit, due to the trivial dependence 
on the interchange of the spin-up and spin-down DOS of 
the lead at lower voltage. It is interesting to notice that 
the TMR reaches a maximum of 18% for asymmetric cou- 
plings. The maximum would increase further as the cou- 
plings become more asymmetric. The TMR ratio given 
by our simple model is consistent with the estimation for 
a resonant structure with Fe electrodes 5 ^ and that of the 
Coulomb-blockade-free double junction model. 54 In mag- 
netic resonant structures, TMR depends not only on the 
DOS of two electrodes as in the single junction case, but 
also on the spectral density of the central well associated 
with G r c c , and thus exhibits complicated dependence on 
the bias voltage. As for the coupling dependence of peaks 
and valleys in the TMR curve, it is associated with the 
sensitivity to the distortion in the spin-up and spin-down 
DOS of the leads. Such a sensitivity strongly depends 
on the coupling strengths and which type of the Fermi 
sea the well levels fall into. In general, electrons in a 
one-band Fermi sea in the weak coupling case can detect 
much better the change in the DOS of the other ferro- 
magnetic lead, and so the TMR develops a peak at the 
boundary between the two distinct types of Fermi sea. 

To summarize this subsection, we have studied the I- 
V characteristics and TMR behavior in a double-barrier 
magnetic structure. It is found that both a peak and 
a shoulder emerge within the resonant region, manifest- 
ing directly the DOS profile of the ferromagnets. This 
finding may provide a new way to measure the degree of 
spin polarization of a ferromagnet. The TMR of reso- 
nant structures exhibits complex dependence on the bias 
voltage. It is either enhanced or suppressed, depending 
on the strengths and symmetry of the elastic couplings 
of the central region to the magnetic leads. 



B. F/I/N/I/S structures 

At an N/S interface a dissipative current in the nor- 
mal metal can be converted into a dissipationless super- 
current in the superconductor via the Andreev reflection 
processi^SiSi Owing to the spin imbalance in the ferro- 
magnet, the Andreev current is suppressed in a F/S 
contact. 36 Blonder, Tinkham and Klapwijk presented 
a one-dimensional model based on the Bogoliubov-de- 
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FIG. 4: TMR versus bias voltage of a magnetic double-barrier 
structure for different couplings at temperature fc_gT = 0.11F. 
The thick lines are the results for TMR, and the thin solid 
and dashed lines represent the tunneling current with paral- 
lel and antiparallel magnetizations, respectively. The other 
parameters are the same as in Fig. 3. 




FIG. 5: A schematic potential profile for a biased DBTS con- 
nected to a ferromagnetic and a superconducting leads. The 
magnetization of the ferromagnetic lead is he and the energy 
gap of the superconductor is |A-r.|. The hatched regions rep- 
resent occupied electron states. A typical Andreev reflection 
process is shown: a spin-up electron above the chemical po- 
tential of the superconductor is reflected as a spin-down hole 
below the chemical potential at the NS interface, and finally 
into the ferromagnetic lead. 



Gennes equation to analyze the transport processes at 
an N / S interface in terms of normal electron transmis- 
sion and Andreev reflection probability, known as the 
BTK theory^ Cuevas et al. in 1996 also uncovered some 
kinds of electron tunneling processes in the N/S quan- 
tum point contacts within the Keldysh NEGF formalism 
starting from a microscopic Hamiltonian^ The scatter- 
ing matrix theory 56 ! 57 and Keldysh NEGF formalism 14 ! 15 
of electronic transport in N — QD — S systems were 
also presented. Quite recently Zhu et alii investigated 
a 2F - QD - S structure using the Keldysh NEGF 
method, obtaining some interesting results. However, in 
the Keldysh NEGF treatment to the N - QD - Si^i 5 . 
or IF — QD — Sf& they always made some assumptions 
that the ferromagnetic magnetization is along the z axis 
and the superconductor order parameter is a real quan- 
tity, and take the wide-band limit. Here we use the cur- 
rent formula (|34|l for F/I/N/I/S systems to investigate 



the resonant Andreev current and I-V characteristics of 
a genuine non-interacting hybrid structure (see Fig. 5) 
beyond the wide-band limit. 

Following similar procedures as in the last subsection, 
we derive the various kinds of full Green functions of the 
normal region for a noninteracting F/I/N/I/S resonant 
structure as 

G^/<(M') = J |;e- te <*-f (47) 
where (X = G r < a /<, g r > a /< and £ r > a /<) 
X(t,0 = P(nnet 



2 ,R/(^)X(t,f) 
R/t ( %) P t(^ c t' 



and 

C# c 8 (e) 



Z'r 1 {e)±-T^{e)± 
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l -T^(e T eV)\ , (48) 
G< c (e) = Gl c (e)[T c f(eTeV)tc(eTeV) 

Substituting the above Green functions into Eq. (|34() . 
we obtain 

T — T A J_ T N 

-Lfns — -Lfns ~t~ 

= ^ J de[f c (e - eV) - f c (s + eV)]T f A ns (s) 
2e 



where 



Ti* 

frts 



+ T y de[/ £ ( £ -ey)-/ TC (e)]T^ s ( £ ), (49) 

iE( r£/ (^ey)G^ c ( £ )) ii+i 

(T^(E T eV)G a CtC (e)) 
V : / i+ii 



i=l,3 



rf f (e-eV)T^(e + eV)\G r c ^ + 

ir 

| ^ -r y|w CjC;34 



rf'( e - e v)rf( e + e v)|& ci34 | a , (50) 



rrN 
1 fns 



\ (r^(£Te7)G^( £ )rf( £ )G» c ( £ )) M 

i=l,3 " 

= i P ^(s)r^{rf ( e - e y)[|G; c;11 | 2 + |g; c;21 | 2 

_2^ Re {G; c;11 Q c;21 }] + rf f (e - eV) 
[l^c.c^sl 2 + |Gc jC;43 | 2 + 

2^Re{G; c . 33 G« c;43 }]}. (51) 

The elements of the Green function matrix & r J^. 111 



LT c,c;12' Ur c,c;21' °c,c;33' °c,c;34 clnu ^ c,cAZ 

from Eq. J3BJ| as 



are derived 



[rf( £ + ey)+r Ks e Ks ( £ )]} 



. (52) 



[if ( e -eV) 



+r*V0O]}{(£ 



2 L i 
1 



+ 



[rf'ce + evo + r 71 "^ (e)]} + 



r^^ s ( £ )i 



»7J 



(53) 



Compared to the work for the N — QD — S system from 
the similar Keldysh formalism^ the derivation of the fi- 
nal current formula (|49[) . (|5UI) and (|51l) from the formal- 
ism we developed is more direct, simple and systematic. 
What we need to do is just some simple matrix algebra, 
while complicated mathematical techniques are needed 
in the derivation of the Green functions in the formalism 
of Sun et al.. 1 — Also the current formula permits us to 
investigate the I-V characteristics within a much wider 
bias voltage region in F/I/N/I/S DBTSs. 

The ferromagnetism and superconductor proximity are 
manifested in the dependence on the magnetization he 
and the magnitude of the superconducting order pa- 
rameter |A-rJ of the full Green functions through self- 
energy matrices. From expressions (|49I5UI51|1 . one ob- 
serves that the current through a F/I/N/I/S resonant 
structure results from different contributions. 13 ' 55 
is the Andreev reflection current: a spin- up/down elec- 
tron/hole associated with spectral weight I\ (e — eV)/ 
r, (s — eV) incident from the ferromagnetic lead is re- 
flected as a spin-down/up hole/electron with spectral 
weight if (e + eV)/ if (e + eV) backward into the 
original lead, and at the same time two electrons in 
the normal region are removed into the superconductor 
as an electron pair with probability \G r c c . 12 \ 2 /\G r c c . M \ 2 . 



Gr 
c,c;ll 

G" r 



G fsl (e) 



c ' c;21 r (r\ 

Gfsi{£) 

E J+F^r 1 + I [rf f (e + eV) + T^g^(e)] 



G 



c.c;33 



Lz c,c;34 — ^c^S 



in which 
Gfsi{e) = {(£ 



(J2 £ +e' 1 +io+ ) 1 1 \Fl ( e F ns (e)] Xf comes from three kinds of physical processes. The 

first and fourth terms in l|51|l represent the contribution 
from normal electron transmission, a spin-up/down elec- 
tron/hole tunnels into the superconductor with proba- 
bility \G r c c -n\ 2 /\G r c c . 33 | 2 ; the second and fifth terms de- 
scribe the 'branch-crossing' process in the BTK theory^ 
a spin-up/down electron/hole in the ferromagnetic lead 
is converted into a spin-down/up hole/electron in the 
superconductor side, with particle pairs of opposite 
spins created in the normal region \G r c c . 21 \ 2 /\G r c C . A3 \ 2 . 
The terms left correspond to the net transfer of elec- 
trons/holes, along with the creation/annihilation of par- 
ticle pairs inside the well and the annihilation/creation of 
pairs into the superconductor lead with probability pro- 
portional to Re{G^ ;11 Q c;21 }/Re{G^ c . 33 G^ ;43 }. At 
absolute zero temperature, the only contribution to the 
current \slf ns for eV < \An\, since in this case p Us {e) in 
Tf^ s becomes zero and then I^ ns = 0. When eV > \An\ 



GfMe) 
G f s2(e) ' 



■ 1 + 4rrf / ' 
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all processes contribute to the current. If one sets he = 
and assumes the wide-band approximation, the current 
formula will reduce to the result obtained by Sun et ali£ 
in the N — QD — S case. 

Assuming a single active level £o in the well, we get the 
following linear-response conductance of the F/I/N/I/S 
system 



Gfns (£o) 
4e 2 



rf / rf / (r w -) 2 /4 



rf'rf' 



(T 



Hs\2 . 



^(r 



T c f? 



(54) 



For the completely polarized ferromagnetic lead, i.e., 
hc/W = 1, and r, = 0, the linear conductance turns 
out to be zero, since there is no state available for the 
Andreev reflected spin-down holes. If the magnetization 
he is zero, = r, ' = L £ (0), the ferromagnetic lead 
becomes a normal metal, and the conductance is reduced 
to 



Gnns\ s o) 



4e^ 
~h 



2T c (Q)T Us 



As 2 



[r £ (o)] 2 + {v ns f 



(55) 



which is the same as the result obtained by Beenakker 
from the scattering matrix approach^ In contrast to a 
single F/S junction^ the conductance of an N/I/N/I/S 
resonant structure is always not less than that of the 
F/I/N/I/S structure, regardless of the value of the 
magnetization he- At £o = the conductance in the 
N/I/N/I/S structure is maximal for symmetric cou- 
plings r £ (0) = L 7 ^, equaling to Ae 2 /h twice that in 
the N/I/N/I/N case. Moreover the line shapes of 
the linear-response conductances (|54() and (|55|l which 
decay as e^~ 4 are not of the simple Lorentzian form 

r c r n /[e 2 + (T c + r n ) 2 /4\. 

Let us analyze further the spin polarization depen- 
dence of Gjns at resonance £o = for different couplings. 
Setting T 7 ^ = Ar £ ^(0), Eq. fg£ evolves into 



Qfr, 



4e 2 4kA 2 
~h~ (k + A 2 ) 2 ' 



(56) 



where k = v/l ~ h 2 c /W 2 is a quantity characterizing 
the degree of spin polarization of the ferromagnetic lead: 
k = 1 for normal metals and k = for completely po- 
larized ferromagnets. When A > 1, the conductance in- 
creases with increasing n, implying that the conductance 
decreases when the degree of spin polarization rises. If 
A < 1, the conductance first increases with increasing 
the spin polarization and then decreases rapidly after it 
reaches its maximum value 4e 2 /h. The critical value of 
k is given by k — A 2 . This interesting result is also ob- 
tained by Zhu et al.^i shown in Fig. 2 of their paper. 

Next we explore the dependence of the Andreev cur- 
rent spectrum on the degree of spin polarization of the 
ferromagnetic lead. From the schematic view of the res- 
onant Andreev reflection processes in F/I/N/I/S struc- 
tures, one immediately becomes aware that the resonant 
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FIG. 6: The Andreev current spectra at zero temperature and 
fixed bias voltage for different spin polarizations (a) he /W = 
, (b) hc/W = 0.75 and (c) hc/W = 0.9, where W = 
2 1 A?j | . The full lines correspond to the results when eV > 
and the dashed when eV < 0. Labels S, D and T are 
used to denote the current peaks arising from the resonant 
Andreev reflections involving single, double or triple levels, 
respectively. Here we assume 20 levels with identical level 
separation 0.5 1 Arc | , and the first level aligns with the chemical 
potential of the superconductor lead when V g = 0. The other 



parameters are T Lt (0) = T Ks = 0.01[A W | 



Andreev current is determined mainly by the applied 
bias voltage V, the ratio of the strength of ferromag- 
netic magnetization to the bandwidth hc/W, and the 
level separation 5e n . We choose a special level separa- 
tion Ae„ = 0.501|Ak|, a case in which at most three 
levels are allowed to fall into the energy gap of the su- 
perconductor lead. For simplicity we assume identical 
level separations and do not consider the influence of the 
bias voltage on the level shift for convenience of compar- 
ison. At fixed bias voltage smaller than the energy gap 
|Atj|, resonant Andreev reflection takes place whenever 
the chemical potential of the superconductor lies just in 
between two levels and there are states available for the 
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reflected electrons/holes. The energy levels e„ — eV g can 
be shifted up and down by tuning continuously the gate 
voltage V g . Therefore, one can expects a series of peaks 
in the Andreev current as a function of V g . 

In Fig. 6 we present numerical results of the Andreev 
current as a function of the gate voltage eV g for differ- 
ent spin polarizations hc/W and different bias voltages 
eV. The cases of positive and negative bias voltage are 
considered to compare the Andreev current contributed 
from electron and hole transmission. Let us first in- 
spect the Andreev current spectra in the N/I/N/I/S 
case (Fig. 6(al)-(6a3)). 15 At a small positive bias volt- 
age eV = 0.1|A-R,|(Fig. 6(al)), a series of peaks labelled 
by S with the same separation as the level spacing is 
observed. These peaks come from the resonant Andreev 
reflection processes by electron tunneling through a single 
level aligned with the chemical potential of the supercon- 
ductor. When the bias voltage is small, there is no pos- 
sibility for two levels to satisfy the resonant Andreev re- 
flection condition. As the bias increases to eV — 0.3\An\ 
(Fig. 6(a2)), there is possibility for two neighboring lev- 
els to lie equally above and below the chemical potential 
of the superconductor. The condition of the resonant 
Andreev reflection involving two levels can be satisfied 
and two-level Andreev reflection also contributes to the 
Andreev current. As a result another series of resonant 
Andreev current peaks labelled by D is observed neigh- 
boring to the original series from single level contribu- 
tions. These D peaks, with the same spacing as the S 
series and 0.25 1 A-r\ away from it, stand out for their dou- 
ble height compared to the S ones. This is because for 
the D peaks two neighboring levels are involved in the 
corresponding resonant Andreev reflection processes, so 
the probability is doubled as compared to the single level 
situation. At a still higher bias voltage eV = 0.6|A^| 
(Fig. 6(a3)) the Andreev current spectrum can be un- 
derstood similarly, the 5-type peaks are now replaced 
by T ones with tripled amplitudes, resulting from three- 
level contributions. Since at eV — 0.6|A-r,|, two addi- 
tional levels near the middle one in alignment with the 
chemical potential also contribute to the Andreev cur- 
rent, so there are three neighboring levels taking part in 
resonant Andreev reflections, making the height of the 
T series three times as that of the S ones. When the 
bias is reversed, the Andreev current becomes negative, 
implying that the Andreev reflection is induced by hole 
transmission. However, the Andreev current spectra re- 
mains unchanged. Since in normal metals the DOS is 
spin degenerate, complete resonant Andreev reflections 
are guaranteed for electrons as well as their hole counter- 
parts. Therefore, except for the sign the Andreev current 
spectra are the same for electron (V > 0) and hole trans- 
mission (V < 0). The above phenomena can also be 
understood from the intuitive diagrams in Fig. 7, with 
different spin-up and spin-down bands of the ferromag- 
nets replaced by identical ones of the normal metals. 

In a F/I/N/I/S resonant structure the Andreev cur- 
rent depends not only on the position of the quantum well 
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FIG. 7: Schematic views of the resonant Andreev reflection 
processes arising from the electron (a) and hole (b) transmis- 
sion from the ferromagnetic lead. The hatched region rep- 
resents states that are occupied by electrons. The Andreev 
current peaks labelled by S in Fig. 6 is originated from the 
Andreev reflection process involving a single level represented 
by the solid line, D double levels by dashed lines and T triple 
levels by dotted lines. The block arrows stand for the process 
at the edges of the superconducting band involving electron 
tunneling through the level located at e n = |An|, which re- 
sults in sharp peaks in the I-V characteristics. 



levels as in the normal case, but also on whether there 
are available states for the backward reflecting holes. The 
Andreev spectra for electron and hole transmission ap- 
pear to be different if the ferromagnetic lead is of large 
spin polarization. For the completely polarized ferro- 
magnetic lead, no Andreev current is expected due to 
the absence of empty states for the returning holes. We 
therefore choose the spin polarizations hc/W = 0.75 and 
hc/W = 0.9 for our purpose, where W = 2\An\, and the 
results are given in Fig. 6 (bl-b3) and (cl-c3). At a 
low bias voltage eV = 0.1|Atj| there are still available 
states for the reflecting holes and the Andreev current 
exhibits the same resonant spectrum as in the normal 
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case, with slightly suppressed peak amplitude. If the bias 
voltage is eV — 0.3|A K |, the spin-down band of polariza- 
tion hc/W = 0.9 moves above the chemical potential of 
the superconducting lead, leaving only the possibility for 
a spin-up electron to be transmitted through a level be- 
low the chemical potential and then reflected backwards 
through the neighboring level above the chemical poten- 
tial to the spin-down band. Hence we can only observe 
the D-type peaks with amplitude half that in the normal 
case in the spectrum (Fig. 6(c2)). As hc/W = 0.75, the 
-D-type Andreev current peaks with half the amplitude 
of the normal case are originated from the same reso- 
nant Andreev reflection processes by spin-up electrons 
going through the states below the chemical potential as 
in the case hc/W = 0.9. This is the reason why we ob- 
serve the amplitudes of the S and TJ-type peaks to be 
nearly the same (Fig. 6(b2)). As the bias voltage eV 
is further increased to 0.6|A-r.|, the spin-down bands for 
hc/W = 0.75 and hc/W = 0.9 shift above the chemical 
potential of the superconductor, and one can no longer 
observe the series of -S-type peaks from the contribution 
of the levels just at the chemical potential. For ferro- 
magnets with high spin polarization, the current arises 
only from the resonant Andreev reflections by the spin- 
up electrons tunneling through the levels 0.5|Ak.| below 
the chemical potential. Thus the Andreev current spec- 
trum ((Fig. 6(c3)) only consists of a series of peaks at 
the positions of the T-type peaks in the normal case. 
Whereas in the small polarization case ((Fig. 6(b3)), the 
Andreev reflections involving two neighboring levels and 
two of three levels contributing to the T-type peaks in the 
normal case can happen, and a series of resonant peaks 
with equal separation 0.25|A-r.| is observed. 

Schematic views of the above resonant Andreev reflec- 
tion processes for electron transmission are given in Fig. 
7(a). When the bias voltage is reversed, the current is 
contributed from hole transmissions, and the situation is 
now very similar to the normal case. The only difference 
is that the amplitudes of the peaks in the ferromagnetic 
case are suppressed due to the reduced DOS for the re- 
flected electrons. The relation among the S-, D- and 
T-type current amplitudes still hold, for which a heuris- 
tic physical picture is given in the Fig. 7(b). When 
the level spacings are not identical, more interesting and 
complicated resonant Andreev current patterns can be 
expected. 15 . However, we can still appreciate and analyze 
them from the intuitive pictures of Fig. 7 whatever the 
Andreev current spectra may be. 

The I-V characteristics of this kind of generic hybrid 
structure is also interesting. It is known that the reso- 
nant Andreev reflection process also contributes to the 
current when the bias voltage eV is greater than |A-r|, 
the energy gap of the superconducting lead. For this 
reason we consider cases in which the level spacing Ae n 
can be either smaller or greater than the energy gap 
| Arc |, as well as when the first level Eq is either below 
or above the chemical potential of the superconductor 
when V = 0. Fig. (8a) and (8b) give the I-V curves for 
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FIG. 8: I-V characteristics of a F/I/N/I /S resonant struc- 
ture at temperature 0.1|A-r| for small level spacing Ae n = 
(a) 0.33|Afj.| below and (b) 0.25|A-r.| above the chemical po- 
tential of the superconducting lead when V = 0. The solid 
curves are the results for the N/I/N/I/S case with symmet- 
ric coupling F Cf (0) = r ns = 0.01|A K |, and the other curves 
for the F/I/N/I/S structure with hc/W = 0.5. 



small level spacing Ae n — 0.5|A-r,| when the first level 
to = -0.33|A K | lies below (Fig. 8a) and e = 0.25|A^| 
above (Fig. 8b) the chemical potential of the supercon- 
ductor lead. As usual we approximate the variation of 
the energy levels with the bias voltage by e„ + 0.5eV. 
There is no substantial difference between the I-V char- 
acteristics of the F/I/N/I/S and N/I/N/I/S systems 
when the level separation is small, as shown by the full 
and dashed lines in Fig. 8. However, the current as a 
function of the bias voltage is strongly dependent on the 
symmetry between the couplings and the level configu- 
ration, especially when the applied bias is positive. The 
peaks in the I-V curves are originated from the resonant 
Andreev reflections, which emerge at some specific bias 
voltages when the resonant Andreev reflection condition 
is satisfied. The irregular current plateaus come from the 
normal particle transmission. There are two types of cur- 
rent plateaus with different heights. Those with higher 
height are determined by particle tunneling through the 
level aligned with one of the edges of the superconducting 



19 



o.o 



c 

Z5 

■H — ' 

Z5 

o 



-0.1 - 



-0.2 
0.1 

0.0 

-0.1 

-0.2 




- NORMAL CASE 

— r L, (0)=r Rs =o.oi/|A R | 
r L, (0)=ior Rs =o.oi/|A R 




NORMAL CASE 

-r Lf (0)=r Rs =o.oi/|A R 



r L, (0)=ior Rs =o.oi/|A R | 



r =ior u (0)=o.oi/|A B 

i' i i 



-10 -8 -6 -4 -2 2 4 6 8 10 

Bias voltage eV/|A | 



FIG. 9: Same as Fig. 8, but for level spacing Ae n = 
2.5|A TC | > W = 2|A K ] when the first level e lies (a) 0.33|A TC | 
below and (b) 0.25|A-7j above the chemical potential of the 
superconducting lead when V = 0. 



energy gap at which the DOS is divergent, represented by 
a block arrow in Fig. 7, while the others are contributed 
from the levels away from the edges. Such an analysis 
can be confirmed by the following simple estimation: If 
e a = -0.33|Att|, then at bias -0.33|A TC | + 0.5eV = 
and 0.17|A-rJ + 0.5eV — 0, the resonant Andreev re- 
flection condition is satisfied and current peaks emerge 
at eV = 0.66|A K | and eV = -0.34|ArJ. For the first 
plateau to appear, we need e n + 0.5eV = ±|Atj| and 
eV > ±\A n \, i.e., eV = 1.66|A TC | and eV = -1.34|A TC |. 
The numerical results are consistent with this simple ar- 
gument (Fig. 8(a)). 

In the case of negative bias eV < — \An\, the levels are 
pushed down gradually inside the energy gap and the 
resonant Andreev reflection condition is eventually satis- 
fied at some bias voltages. Since both the normal particle 
tunneling and the resonant Andreev reflection contribute 
to the current, one observes a series of equally spaced res- 
onant peaks superimposed onto the plateaus with identi- 
cal widths. When the couplings become asymmetric (one 
of the couplings is weakened), both the Andreev current 
peaks and current plateaus are suppressed. This is be- 
cause the normal and Andreev current is proportional to 



the product of the couplings r £ ^(0)r Ks , as can be found 
in the formulae (|49|l . (|50|l and (|51|l . It has already been 
found analytically by Sun et a! 15 that the Andreev reflec- 
tion probability is maximized for symmetric couplings in 
N — QD — S structures and decreases rapidly with the 
increasing coupling asymmetry. In the special asymmet- 
ric coupling case where T ns << r £ ^(0), the sharp DOS 
at the edges of the superconductor energy gap can be 
discerned clearly whenever particles transmit through a 
level aligned to the edge (see the dotted lines in Fig. 
8) . This can be easily understood from the Breit-Wigner 
formula for resonant normal electron transmission at ves- 
ona,nceT = 4r c f(0)r ns p c fp ns /{r c f(0)p c f+r Ks p ns } 2 . 
The transmission probability T depends on the ratio be- 
tween the level-widths T cf (0)p cf and T Ks p ns . At the 
edges of the gap, p Ks is divergent and thus the coupling 
constant should be small enough to balance these two 
level-widths to guarantee high transmission. Notice that 
the ferromagnetic feature can not be observed in the I- 
V curves of F/I/N/I/S structures when the level spac- 
ing is small. It can only be displayed when the level 
separation is greater than the band-width of the ferro- 
magnets. In Fig. 9 we present the results of I-V re- 
lation of this case for different couplings and different 
level arrangements. As expected, the degree of spin po- 
larization of the ferromagnetic lead is reflected in the 
I-V curves. In the normal case and when eV > 0, 
the current first develops a resonant Andreev peak at 
— 0.33|A7j|+0.5eV = and exhibits usual resonant peaks 
in the double-barrier structure after a narrow peak with 
width determined from — 0.33| A75.I + Q.^eV = \A-r\ and 
-0.33|A TC | + 0.5eV = eV - W. When the bias is re- 
versed, the current displays a plateau from the position 
— 0.33| A75.I + 0.5eV = — |A^|, superimposed by some 
weak Andreev peaks arising from the resonant Andreev 
reflections for |e| > \A n \. If e = 0.25|A TC | when V = 0, 
one observes an I-V characteristics similar to that of a 
usual magnetic DBTS in the case of positive bias V > 0, 
and in the negative bias case a resonant Andreev peak 
at 0.25 + O.beV — and a plateau-peak structure simi- 
lar to the case e = -0.33|A K |. For F/I/N/I/S struc- 
tures, the current develops both a resonant peak as well 
as a shoulder which are more prominent in the high bias 
case for eV > 0, while they keep nearly the same as 
in the N/I/N/I/S structure when the bias is negative. 
This picture is violated when the couplings become asym- 
metric. As in the case of small level spacing, the sharp 
DOS at the edges of the energy gap is also prominent 
in the I-V characteristics when the coupling to the su- 
perconductor side r Ks is much smaller than that to the 
ferromagnetic side r -C ^(0). The spin polarization of the 
ferromagnet can thus be measured when the couplings 
to the ferromagnetic and superconducting leads are sym- 
metric, as suggested from the comparison of the dashed, 
dotted and dash-dotted lines in Fig. 9. It is noted that 
the I-V characteristics of F/I/N/I/S resonant structures 
are qualitatively right, because the idea about the band- 
width of the superconductor is somewhat vague within 
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its semiconductor models 

In summary, we have investigated in this subsec- 
tion the Andreev current spectra and I-V relations of 
a F / 1 /N / 1 / S resonant structure in detail. Interesting 
dependence on the ferromagnetic spin polarization of the 
linear conductance is discussed in terms of analytic ex- 
pressions, given by Eqs. (|54|l and l|56|l . Our results 
demonstrate that the peak structure of the Andreev cur- 
rent as a function of the gate voltage is determined by the 
applied bias (both the value and sign) and the degree of 
spin polarization, which differs substantially from the re- 
sults under the wide-band approximation. The I-V char- 
acteristics, closely associated with the level arrangement, 
coupling symmetry, bias sign and spin polarization, can 
be employed to characterize the density of states (DOS) 
of both ferromagnets and superconductors by tuning the 
coupling strengths. 



C. S/I/N/I/S structures 

The discovery of the Josephson effected has provoked a 
lasting research interest in the properties of the DC-and 
AC- Josephson current in mesoscopic S/N/S junctions in 
thirty years^LSi When the width of the normal region 
is smaller than the coherence length, electron pairs can 
coherently tunnel from one superconductor to the other, 
inducing a phase-dependent DC current even when the 
bias is zero. Early in 1963, Ambergaokar and Baratoff 61 
derived a useful formula for the supercurrcnt in S/I/S 
junctions with the help of the Gor'kov Green functions. 
In the 1990's investigations on the mesoscopic S/N/S 
junctions became timely^ due to the advances in exper- 
imental techniques. In most of these work, the scatter- 
ing matrix method based on the Bogoliuboiv-de Gcnncs 
(BdG) equation^ is commonly used. In a S/N/S junc- 
tion, the Andreev reflections at the N / S boundaries con- 
fined the quasi-particle inside the normal region, result- 
ing in the bound states sensitive to the superconducting 
phase difference of the two superconductors. 32 Impurities 
inside the normal region altering the quasi-particle wave 
interference^ and the asymmetry between the two en- 
ergy gaps can also modify the Josephson current^ Glaz- 
man and Matreev^Si and Ishizaka et alS have studied 
the influence of the Coulomb interactions on the Joseph- 
son current in S — QD — S systems. Research were also 
conducted on the DC Josephson current in noninteract- 
ing symmetric S — QD — S structures by Beenakker us- 
ing the scattering matrix approach;--? and by Lin's Group 
from the Keldysh NEGF method. ?£ However, these inves- 
tigations are restricted to the symmetric case-the same 
couplings and energy gaps. Motivated by this limitation, 
we investigate in this subsection the DC Josephson cur- 
rent through a general S/I/N/I/S resonant structure, 
in order to reveal the dependence of the Josephson cur- 
rent on the energy gaps. Results for the AC Josephson 
current will be reported elsewhere. 

Neglecting the interaction effects in the central region 




FIG. 10: A schematic potential profile for a biased DBTS 
attached to two superconducting leads with order parameters 
Ac and An- Here the hatched region represent occupied 
electron states. 



we obtain an unexpectedly simple form of the DC Joseph- 
son current formula for a general S/I/N/I/S system 
fromEq. © 
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where we assume in general \Ac\ < |Atj|, and e p , the 
energies of the discrete Andreev bound states, are the 
poles of the spectral function G r s (e) 
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It is seen from Eq. H57JI that the DC Josephson cur- 
rent T sns in a general S/I/N/I /S system has contribu- 
tions from three different scattering processes: T\ results 
from the resonant Josephson tunneling through the dis- 
crete Andreev bound states given by G r s (e) — within 
|e| < |A,c|; X 2 from the quasi-particle escaping through 
broadened levels from the normal region to the weaker su- 
perconductor side, i.e., |A,c| < \e\ < \A-ji\; and T3 from 
quasi-particle tunneling from the normal region to both 
superconductors |e| > |Atj|.S One can consider Eq. I|57|) 
as an extension to the asymmetric case of the Beenakker's 
result for symmetric S — QD — S systems In ad- 
dition one can check after simple algebra that the DC 
Josephson current formula i|57|) for the general asymmet- 
ric S/I/N/I/S system reduces to the known result for 
the symmetric casei^SiSLS 

Now we consider the simplest situation in which there 
is only a single active level £0 in the central normal re- 
gion. It is expected that the resonant Josephson scatter- 
ing via the Andreev bound states e p will dominate the 
DC Josephson current. In Fig. 11 we plot this quantity 
at zero temperature, calculated from Eq. I|57(l as a func- 
tion of the single level energy £ f° r symmetric couplings 
T Cs = T 718 = 0.01|A/;|. However, the superconducting 
energy gaps are allowed to be asymmetric. The supercon- 
ducting phase difference is chosen &sip s = ir/2. The total 
current labelled by / in Fig. 11(a) has a resonant peak 
when the single level is aligned with the chemical poten- 
tial of the superconductor, i.e., £0 = 0, resulting from the 
constructive interference between the forward Andreev 



state +£p and backward Andreev state 



-e p . Inspecting 



the current components for three different scattering re- 
gions \e\ < |A £ |, |A £ | < \e\ < |A W | and \e\ > Ar 1 ^ 
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FIG. 11: The DC Josephson current vs £0 for (a) different 
scattering processes, and (b) different energy gaps, where we 
take symmetric couplings as r £s = T TCs = 0.01|A £ |. The 
superconducting phase difference is ip s = 7r/2. 



labelled respectively by 1, 2 and 3, one finds that the cur- 
rent component 1 contributed by discrete Andreev levels 
makes the major contribution to the DC Josephson cur- 
rent. It possesses one peak at £0 = plus two side peaks 
pinned at £0 = ±|A,c| which are offset by two peaks 
in the current component I2, in which two additional 
wider side peaks cancelled exactly by X 3 . The side peaks 
come from the abnormal superconductor DOS singular- 
ities at the edges of the energy gap(s), and the exact 
cancellation of these side peaks in the total Josephson 
current is due to the fact that both electron-like and hole- 
like excitations can escape through the active level from 
the superconductor to the central normal region. This 
is equivalent, mathematically, to the vanishing residues 
of the spectral function G r s ! 67 In contrast to the sym- 
metric case \&c\ = \^n\, the DC Josephson current is 
slightly enhanced in the asymmetric case \Ac\ < |A-r.|, 
as demonstrated in Fig. 11(b). This result differs from 
the usual S/N/S structure, where the current is greatly 
enhanced^ The reason is that the energy levels lie within 
|e| < |A,c|, and then resonant Josephson tunneling dom- 
inates the DC current, as can be seen more clearly in 
the current-phase relation in Fig. 12. The significant 
enhancement of the DC Josephson current may be ob- 
served in a resonant S/I/N/I/S system with many lev- 
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FIG. 12: The current-phase relation for a resonant 
S/I/N/I/S structure, where the insets are enlargements of 
the current component Z3 . Here the total DC Josephson cur- 
rent T sns and its components X\, X2 and X3 are represented 
by the thick solid, thin solid, dashed and dotted lines, respec- 
tively. 



els located inside the region |e| < |A,c|, since in this case 
so many levels are active to contribute to the Joseph- 
son current. Fig. 12 shows the resonant DC Josephson 
current at zero temperature as a function of the super- 
conducting phase difference (p 8 for general asymmetric 
S/I/N/I/S systems. When e = 0, i.e. the single level 
is exactly located at the position of the superconducting 
chemical potential, the DC Josephson current I sns (ip s ) 
and its component X\ vs the superconducting phase dif- 
ference (p s is of the sawtooth shape (Fig. 12 (a) and 
(b)), no matter how big the difference between the two 
superconducting energy gaps is. The physical origin of 
the sharp discontinuity at (p s = ir is the same as in the 
usual S/N/S junctioniSiSi the Andreev levels e v and 

— e p determined from the spectral function G r s {e) = 
at (f s — interchange their position at ip s — n, produc- 
ing a discontinuity in the DC Josephson current. Unlike 
the usual S/N/S structure, the supercurrent-phase rela- 
tion has a weak dependence on the asymmetry between 
the superconducting energy gaps in resonant S/I/N/I/S 
structures. This is because the main contribution to the 



DC Josephson current in resonant structures with a sin- 
gle level is mainly from the Andreev refection processes 
inside the region |e| < \Ac\, where the energy of the An- 
dreev bound states has a trivial dependence on the energy 
gap |Atj|, as can be seen from the spectral function as- 
sociated with Gg(s). When the single level moves away 
from the position of the chemical potential £0 = 0.2|A,c| 
(Fig. 12(c)), or the coupling to the superconductors be- 
comes asymmetric T Ks — 4T Cs (Fig. 12(d)), the abrupt 
jump at ip s — 7T is smoothed out and the current-phase 
approaches the sinusoidal relation. The transmission 
probability becomes smaller when the level moves away 
from the chemical potential or the elastic couplings be- 
comes asymmetric. As a result the link between these two 
superconductors becomes wicker, and thus the current- 
phase relation Z(tp s ) oc smip s is expected^ As the sin- 
gle level moves away from the chemical potential or the 
couplings become unequal, the DC Josephson current is 
significantly suppressed, with the component X\ from the 
discrete spectra suppressed while the component I2 from 
the continuum spectra \^c\ < | £ | < |A^| enhanced. The 
suppression of X\ is due to the decrease in the resonant 
Josephson tunneling probability which is originated from 
the violation of the constructive interference between the 
wave-functions of the Andreev levels, while the enhance- 
ment of I2 results from the fact that the Andreev levels 
are pushed towards the region |A^| < |e| < |A-r| and 
thus the leaky probability is increased. 

In summary, we have shown in this subsection that 
the DC Josephson current in an asymmetric-gap reso- 
nant S/I/N/I/S structure with a single level is slightly 
enhanced in contrast to the symmetric-gap case. The 
current-phase relation is closely related to the position of 
the single level and the symmetry between the couplings 
to the two superconducting leads. 



IV. CONCLUDING REMARKS 

We have developed a unified theory of electronic trans- 
port in a general two-terminal hybrid nanosystem, in 
which each lead can be either a ferromagnet or a super- 
conductor. Within the Keldysh NEGF formalism, the 
current is expressed in terms of the local properties of 
the central interacting region, 

Qr,a,</> and the equilib- 
rium distribution functions of the leads f 7 . The ferro- 
magnetism and superconducting proximity are treated 
on the same footing, incorporated into the tunneling 
Hamiltonian and the self-energy matrices after introduc- 
ing a four-dimensional Nambu-spinor space and perform- 
ing appropriate Bogoliubov transformations. With the 
help of some unitary rotation and phase matrices, one 
can demonstrate analytically the gauge invariance of the 
general current formula H31[). and simplify it to the Meri- 
Wingree forms for specific structures. For some quanti- 
ties, such as the chemical potential, magnetization orien- 
tation and the superconducting order parameter phase, 
only their relative value appears explicitly in the expres- 
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sions of current. Moreover, resonant tunneling, strong 
electron correlations (Coulomb blockade, Kondo reso- 
nance etc.), ferromagnetism and superconductivity prox- 
imity effect can be investigated in a unified transport 
theory without introduction of any ad hoc assumptions. 
In addition, the energy- and bias-voltage-dependence of 
the level-width functions and distribution functions en- 
ters into the current formula in a strict and natural man- 
ner, allowing us to explore the I-V characteristics in the 
large bias limit. However, the disadvantage of applying 
first the Bogoliubov transformation for the ferromagnetic 
lead Hamiltonian is that we can only obtain the expres- 
sions for the sum of the spin-up and spin-down current, 
while the information about the spin components of the 
current is lost. 

Applying the current formulae to the simplest DBTS 
where the interactions are ignored some interesting trans- 
port properties are revealed if one takes into considera- 
tion the finite energy band structure of the ferromagnets. 
We have reported on the current flowing through a non- 
interacting zero-dimensional central region, thus the re- 
sults obtained are qualitatively right for a 2D quantum 
well with the attached emitter and/or the collector be- 
ing ferromagnetic. In addition, we did not consider the 
spin-flip process due to the interfacial scattering or the 
existence of paramagnetic impurities inside the barrier. 
It is known that the spin-flip process may reduce the 
magnitude of the tunnel magnetoresistance of a F/I/F 
junction^ Also, we can expect that the Andreev cur- 
rent spectrum in F/I/N/I/S structures will be modified 
to a great extent in the presence of the spin-flip process, 
since the Andreev reflection may be enhanced with the 
assistance of such process. 

InF-QD-F, F-QD-S or S-QD-S systems, 
electron-electron interactions inside the QD are impor- 
tant and thus one should consider the many-body cor- 
relation effect. One example is the Kondo effect at low 
temperatures^ In such a circumstance, one has to calcu- 
late the full Green functions of the QD in the presence of 
electron-electron interactions, taking into consideration 
the couplings between the QD and the leads. We are 
aware of three recent preprints 7 ^ on the Kondo physics 
in F — QD — F systems in which the wide-band approxi- 
mation is used. 16-72 Such a simplification of the ferromag- 
netic DOS may lead to even spurious results in I— V char- 
acteristics. However, interesting and even unexpected 
Kondo resonances in these systems may arise with the 
full consideration of the finite energy band structure of 
the ferromagnetic leads. 
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APPENDIX A: DERIVATION OF THE 
SELF-ENERGY MATRICES S r '?; </:> 

7/7 s 



In this Appendix we derive various kinds of self-energy 
matrices for the elastic couplings between the central re- 
gion and the ferromagnetic and superconductor leads. 

First, we calculate the retarded/ advanced self-energy 

matrix Tf^ nm (ti,t2) arising from the coupling between 
the central normal metal and the ferromagnetic lead 
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(h,t 2 



= Ev^ (tl)g ;/; 7/t ( tl , i2) v^( i2 ) 



The sum over momentum k can be converted to energy 
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Similarly, after transforming the momentum sum ^ fe 
into an integral J dskpj^ (£k), where pl^ is the normal 
state of the superconductor, we obtain the self-energy 
matrix due to the coupling of the central region to the 
superconducting lead 

H-nmihM) = J2 V lnHh)stk,-y S k(h,t2)VZ(t2) 
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where we have defined the complex level- width matrix as 
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Here we have defined a complex superconducting DOS, 
extending to the forbidden region in the usual BCS the- 
ory | £ | < |A 7 | , inside which the Andreev reflection pro- 
cesses arise, as shown in the Blonder-Tinkham-Klapwijk 
theory^ When |e| < |A 7 |, the quasi-particle density of 
states q is purely imaginary, indicating evanescent states 
in the gap which eventually decay into the pair conden- 
sate. The quasi-particle density of states of the super- 
conducting lead 7 is defined as 
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In deriving Eq. I|A4(1 we have used the following equal- 
ities 
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Note that we have chosen different complex half-planes in 
the contour integrations, in order to guarantee that the 
retarded/advanced self-energy matrices satisfy E 7S (s) = 

The lesser/greater self-energy matrix defined by 

s</> nm ( tl , t2 ) =Evr /st (ii)gS, 7//s (ii,i 2 )vii /8 w 

are obtained in a similar way 
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P(^c*2 + y ), 
where the real level-width matrix is defined as 
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(A10) The Fermi distribution matrix in the Numbu-spinor 
space takes the following form 



|A 7 | 
e— c 



pT s (e-c) p"< s (s + c) J 



f 7 (e=Fc) = 

//(e-c) 
/(e + c) 



V 






f(e-c) 




f 7 (e)=f 7 ( £T 0). 







f(£ + c) 



(All) 
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